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IKTRCDUCTICN 
Scope 

This  report  outlines  the  detailed  structure  of  the  earth's 
dagnetic  field,  indicates  briefly  the  major  phenomena  associated 
with  the  field,  and  considers  some  possible  physical  models  of 
the  field.   I(jiov7ledg:e  of  the  field  structure  and  the  major 
magnetic  phenomena  becomes  increasingly  im.portant  as  man 
ventures  into  the  environment  of  space.   Physical  models,  if 
simulated,  would  be  valuable  in  solving  such  problems  as  the 
origin,  spatial  distribution,  and  future  behavior  of  the 
geomagnetic  field.  Also,  many  of  the  processes  implicit  in  the 
generation  of  the  earth's  field  may  be  valid  in  the  production 
of  magnetic  fields  in  the  sun,  stars,  and  galaxies. 

As  a  starting  point  of  the  discussion,  surface  measurements 
of  the  earth's  field  are  described,  and  spherical  harmonic 
analysis  is  employed  to  obtain  a  magnetic  potential  function. 
From  the  potential  function  the  contribution  of  internal  and 
external  sources  to  the  geomagnetic  field  is  determined.   Next, 
a  description  is  given  of  the  various  short  term  magnetic 
fluctuations  that  are  superimposed  upon  the  earth's  mean  field. 

An  examination  of  the  results  of  the  spherical  harmonic 
analysis  leads  to  several  models  of  the  earth's  field;  these 
models  range  from  a  sim.ple,  gross  representation  to  a  complex, 
accurate  representation.   The  simplest  two  models  describe  the 
earth's  field  in  terms  of  a  centered  dipole  and  an  eccentric 
dipole.  A  third  model  adds  to  the  internal  diDole  field  the 


residual  or  non-dlpole  part  of  the  internal  field.   The 
residual  field  is  noteworthir  since  it  seems  directly  corre- 
lated with  the  secular  or  long  terra  variations  in  the  earth's 
magnetism.   The  field  due  to  external  sources  is  added  to  the 
Internal  field  to  produce  a  final  model;  this  external  field 
appears  to  be  associated  vjith  the  transient  magnetic  fluctu- 
ations. 

Most  theories  advanced  in  the  past  to  explain  the  origin 
of  the  earth's  internal  field  stand  discredited  today.   The 
theory  most  nearly  consistent  with  the  present  knov/ledge  of  the 
earth's  interior  is  the  dynamo  theory.  However,  an  exact  so- 
lution for  a  self-regenerating  internal  dynamo  seems  impossible 
since  so  little  is  knov/n  about  the  boundary  conditions,  and 
since  the  ec_uations  describing  the  process  are  non-linear.  .  Some 
qualitative  information  on  the  djmamo  process  can  be  gained  from 
a  dim.ensional  analysis  of  the  basic  equations  of  magneto -hydro- 
dynamics. 

The  concluding  portion  of  this  report  deals  briefly  v^ith 
some  of  the  phenomena  associated  v;ith  the  external  portion  of 
the  geomagnetic  field.   Possible  theories  on  the  cause  of  mag- 
netic storm.s  are  discussed,  particularly  those  involving  an 
emitted  stream  of  particles  from  the  sun.   Some  of  the  possible 
interactions  betv/een  the  geomagnetic  field  and  these  solar  par- 
ticles are  considered. 

Historical  Background 

The  discovery  in  Asia  Minor  of  a  strange  mineral  having 


the  povfer  of  attracting  ferrous  objects  was  the  subject  of 
great  mystery  and  speculation  to  the  ancient  Greeks.   Bits  of 
this  mineral,  or  magnets  (so  named  after  the  region  Magnesia 
where  deposits  of  magnetic  oxide  of  iron  v;ere  found),  were 
often  exhibited  by  certain  mystic  brotherhoods  as  a  proof  of 
superns.tural  skill.   The  Greeks,  however,  had  no  knowledge  of 
the  poles  of  a  mag;iet,  or  of  the  directive  property  shovm  by 
a  magnet  free  to  move  in  the  earth's  field. 

In  the  eleventh  centur3^  Chinese  \\T7i tings  reported  ships 
using  magnets  mounted  upon  floating  straws.   Some  historians 
believe  that  this  knovrledge  of  the  magnet's  north-seeking 
property  may  have  come  from  Arab  or  Persian  navigators. 

First  mention  of  the  compass  in  European  literature  was 
made  in  1187  V7ith  the  r^ublication  of  a  treatise  by  Alexander 
Neckham  describing  the  action  of  the  m.agnetic  needle  (Mitchell, 
1932) .   The  compass  needle  was  soon  to  be  described  by  Scandi- 
navian sailors  as  the  lodestone  or  "leading  stone".  Magnetic 
dip  or  inclination  was  first  noted  in  1576  when  Robert  Norman, 
a  practical  seaman  and  instrument -maker,  published  a  discourse 
entitled,  "A  newe  discovered  secret  and  subtill  propertie  con- 
cernying  the  decllnying  of  the  Needle,  touched  therewith  under 
the  Plaine  of  the  Horizon" . 

Gilbert  (1600)  postulated  that  the  earth  was  a  magnet,  and 
in  order  to  demonstrate  this  idea  he  shaped  a  sphere  of  lode- 
stone  and  traced  the  field  pattern  using  dip  needles.  Next 
followed  Gellibrand's  announcement  in  1634  that  "variation  is 
accompanied  by  a  variation",  i.e.,  the  angle  the  compass   ,  ' 


needle  made  v;ith  true  north  varied  v.'ith  time.   It  vras  scon 
afterv:ards  reco£nized  hy   Edmund  Kalley  (1705)  that  a  systematic 
survey  of  masnetic  declination  over  the  oceans  v;ould  greatly 
benefit  navigators.   Therefore,  the  first  survey  of  this  t^q^e 
v;as  made  in  the  Atla.ntic  Ocean  during  169S-17C0  by  the  saili-ng 
vessel  Paramour  Fink.   The  data  thus  accumulated  was  published 
as  the  first  isomagnetic  chart. 

By  using  the  potential  theory  developed  by  Laplace  and 
Poisson,  Gauss  (1839)  made  the  first  mathematical  analysis  of 
the  earth's  field.   Since  the  magnetic  survey  data  available 
for  this  analysis  failed  to  cover  extensive  areas  of  the  earth' s 
surface,  G-auss  realized  the  importanc©  of  a  complete  magnetic 
survey.   For  this  purpose  in  1840  a  world-'i-:lde  system  of  mag- 
netic observatories  was  founded. 

Careful  magnetic  measurements  taken  at  these  observatories 
indicated  a  wide  number  of  fluctuations,  both  regular  and  ir- 
regular, in  the  geomagnetic  field.  Particularly  interesting 
v/ere  the  irregular  disturbances  termed  magnetic  storms;  these 
magnetic  storms  were  scon  correlated  with  other  natural  phenom- 
ena such  as:  (1)  earth  currents  f levying  in  a  telegraph  system, 
(2)  the  eleven  year  cycle  in  sun  spot  activity,  and  (3)  auroras 
in  the  polar  regions. 

Further  physical  Investigations  substantiated  the  corre- 
lations just  described.  Birkeland  (1896)  v;as  able  to  produce 
auroras  in  his  laboratory''  by  using  a  magnetized  steel  sphere 
v/ithin  a  vacuumi  discharge  chamber.   The  calculation  of  charged 
particle  paths  in  the  geomagnetic  field  was  performed  by 


Stormer  (1907);  the  results  tended  to  support  the  general  fea- 
tures of  the  auroras.   To  explain  1-Iarconl's  successful  trans- 
mission of  radio  signals  across  the  Atlantic,  Kennelly  s.nd 
Kea.viside  (1902)  proposed  independently  the  existence  of  an 
electrically  conducting  layer  in  the  upper  atmosphere.   This 
layer,  the  ionosphere,  v:as  shown  to  exist,  and  its  behavior  was 
correlated  with  magnetic  storm  activity.  A  theory  proposed  by 
Chapman  and  Ferraro  (1951)  explained  magnetic  storm  activity  in 
terms  of  a  shower  of  emitted  solar  particles  that  tended  to 
distort  the  geomagnetic  field  upon  their  arrival  at  the  earth. 

Repeated  mathematical  analyses  of  the  geomagnetic  field 
have  shovrn  a  general  westward  drift  with  time  of  the  non-dipole 
portion  of  the  field.  Also,  studies  of  the  remanent  magnetism 
in  rocks  suggested  that  the  earth's  magnetic  poles  were  once  at 
positions  far  removed  from  their  location  today,  and  that  pos- 
sibly the  earth's  magnetic  polarity  has  undergone  reversals  in 
past  ages.   These  facts  give  considers-ble  support  to  a  dynamo 
theory  proposed  by  Elsasser  (19-^6)  and  Bullard  (1949a)  as  an 
explanation  of  the  origin  of  the  internal  geomagnetic  field. 

In  recent  years  rocket  and  satellite  measurements  taken 
8.t   distances  out  to  ten  earth  radii  from  the  earth'  s  surface 
have  detected  a  belt  of  charged  particles  trapped  in  the  geo- 
m.agnetic  field.   This  Van  Allen  belt  has  been  artificially 
created  by  the  detonation  of  a  nuclear  device  in  the  earth' s 
high  atmosphere,  as  in  the  Argus  experiment.   Such  a  belt  of 
radiation,  whether  due  to  natural  or  man  made  causes,  is  a 
definite  hazard  to  man's  operation  in  space  and  is  presently 
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the  subject  of  intensive  investigation. 

DESCRIPTICN  OF  THE  EARTH'S  LIAC-ITETIC  FIELD 
Definition  of  Terms 

Two  standard  systems  may  be  used  to  specify  the  magnetic 
field  at  any  point.   The  first  system,  of  primary  use  to  navi- 
gation and  surveying,  defines  the  total  magnetic  intensity 
vector  F  in  terms  of  the  magnetic  elements  K,  D,  and  I.  H  is 
the  magnitude  of  the  horizontal  component  of  F  and  is  consid- 
ered always  positive.   The  vertical  plane  thj?ough  F  is  called 
the  local  magnetic  meridian.   Then,  at  any  point  P,  the  mag- 
netic declination  D  is  the  angle  between  the  magnetic  m^eridian 
and  the  geographical  meridian.  East  declination  is  considered 
positive,  and  v;est  declination  is  negative.   The  angle  F  m.akes 
v;lth  the  horizontal  plane  is  I,  the  magnetic  inclination  or 
dip;  I  is  positive  if  F  is  inclined  dovmwards,  and  negative  if 
F  points  upwards. 

The  second  system,  gives  F  in  terms  of  its  rectangular 
components  X,  Y,  and  Z,  where  X  is  taken  along  the  intersection 
of  the  horizontal  plane  and  the  geographical  meridian;  if 
northward,  X  is  positive,  if  southward,  negative.   Y  is  taken 
in  the  horizontal  plane  transverse  to  the  geographical  m.eridian, 
and  is  positive  if  eastvrard,  and  negative  if  westward.   Z  is.  the 
vertical  component,  and  is  positive  if  dov/nward,  negative  if 
upward . 

From  Fig.  1,  the  following  relations  exist  among  the 
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FIG.  I.    ELEMENTS     OF     THE     EARTHS     MAGNETIC    FIELD 
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na^netic  elements: 

X  =  PI  cos  D 

Y  =  H  sin  D 

Z  =  F  sin  I 

f2  =  x2  +  y2  +  z2  =  h2  +  z2 

H  =  (X2  +  y2)s 

D  =  tan"^  Y/X 

I  =  tan-^  Z/H  (1) 

The  linits  of  magnetic  force  in  geomagnetic  measurements 
are  the  gauss  (f)  and  the  ga.mjna'  (^),  the  gamma  being  defined 
as  10"^  gauss.  Strictly  speaking,  the  unit  of  magnetic  force 
is  the  oersted,  and  the  unit  of  magnetic  flux  densitj'-  is  the 
gauss.   Hov/ever,  the  permeability  of  most  non-ferrous  materials 
is  unity  in  the  cgs  system,  and  little  error  is  introduced  by 
tliis  discrepancy  in  term.inology. 

The  geomagnetic  field  is  not  constant  and  exhibits  two 
distinct  types  of  change,  transient  fluctuations  and  long-term 
secular  changes.   Transient  variations,  irhich  are  discussed  in 
a  later  section  of  this  report,  arise  largeljr  from  causes  ex- 
tern?,l  to  the  earth.   The  secular  variation  extends  over  a 
period  of  many  years  or  even  centuries,  a.nd  is  due  to  cp.uses 
v;ithin  the  earth. 

1-ieasurements  taken  of  the  magnetic  elements  H,  D,  I  or  X, 
Y,  Z  at  a  Fiagnetic  observation  point  are  given  in  terms  of  m.ean 
annual  values  by  a^veraging  out  all  transient  variations.   The 
results  of  a  magnetic  survey  to  determine  the  geomagnetic  field 
over  the  earth's  surface  are  best  illustrated  by  isomagnetic 


charts,  i.e.,  maps  having  lines  drat'/n  tlirouch  all  points  at 
v:iiich  a  given  magnetic  element  has  the  same  value.  A  chart 
giving  lines  of  equal  magnetic  declination  D  is  an  isogenic 
cha.rt,  one  giving  lines  of  equal  inclination  I  is  an  isoclinic 
chart,  and  those  giving  lines  of  equal  X,  Y,  Z,  H,  or  F  are 
isod^Tiamic  charts.   The  annual  secular  variation  of  each  mag- 
netic elem.ent  can  "be  charted  also  in  an  isoporic  chart,  v;here 
points  having  the  same  annual  change  in  a  given  magnetic  ele- 
ment are  connected  hy  lines,  _  ■ 

Spherical  Harmonic  Analj'-sis  ■ 

Within  any  region  of  space  free  of  magnetic  or  current 
distributions,  the  follov/ing  relations  hold  for  the  magnetic 
intensity  P, 

curl  F  =  0 

F  =  -grad  V  (2) 

so  that  V  is  the  magnetic  scalar  potential  satisfying  Laplace's 
equation 

V^V  =0  (3) 

Let  the  set  of  spherical  co-ordinates  (r,  0,  e)  be  defined, 
Inhere  r  is  the  length  of  the  radius  vector,  0  is  the  centrs.l 
azimuthal  angle  or  east  longitude,  and  6  is  the  polar  angle  or 
co-latitude  (Pig.  2). 

For  this  co-ordinate  system  Laplace's  equation  has  two 
forms  of  solution.   If  all  the  magnetic  matter  is  inside  a 
sphere  of  radius  a,  and  none  outside, 
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FIG.   2.       CO-ORDINATE    SYSTEM. 
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1 


Vi  =  a  >  (a/r)^-^^  Sn(e,0)  r  >  a         (4) 


n=0 
If  all  the  magnetic  matter  is  outside  the  sphere,  and  none 

inside. 


: 1 


Ve  =  a  >  (r/a)^  sg(e,0)  r  <  a        (5) 


n=0 
where  Sn(6,0)  is  the  spherical  surface  harmonic  function  of 
degree  n  given  by 

^n  =  6n>0  ^n(©) 

n 

+  /  ^Sn,ni  cos  m0   +  hn,iii  sin  m0)   Pn,m(©)         (6a) 
in=l  (gn  ra  a^^  hn,m  ccnstants) 

Pn(©)  "being  the  Legendre  function  of  degree  n,  and  Pn,in{6)  the  . 
associated  Legendre  function  of  degree  n  and  order  m. 


If 


?n,!n(6)   is  integrated  over  the  region  of 


orthagonality,  then 

1 

/  r^n,in(©)l^  ^(cos  e)  =  2/(2n+l)  (n-nn)!  (7) 

7  (n-in)l 

-1 

Thus  for  any  degree  n  the  mean  square  values  of  Pn  m(6)  have 

very  diverse  orders  of  magnitude.   For  example,  the  ratio  of 

squares  of  P4  1  and  P4  4  is  1  :  2016, 

Schmidt  (1935)  devised  a  set  of  partially  normalized 

Legendre  functions  to  be  used  in  a  spherical  harmonic  analysis 

of  the  geomagnetic  field.   These  functions  ?g(0).are  defined  by 
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Pg(G)    =  Pn(e) 

P8(e)   =     [2(n-m)I/(n+m)I  ]-^  Pn,in(e) 
Sn  may  ncv;  be  re\rrltten 

n 


(8) 


.0^0 


m 


m 


^^l/ 


Sn  =   BnPn(e)    +     >  (Sn  cos  iii0  +  h^  sin  iii0)   Pg(e) 


(6b) 


Ehrl 


SO  the  average,  taken  over  the  surface  of  a  sphere,  of  any  term 

in  Sn  is 

n  2rr 


1/477 1     )  [Pg(e)||jg]ni0]2   g^j^  Q  ^Q^0  ^  i/(2n+l) 


(9) 


0  0 


Therefore,  the  coefficients  g^*  and  h^  in  S^  now  indicate  the 
approximate  orders  of  magnitude  of  the  corresponding  terms. 

In  the  most  general  case  magnetic  material  ma^r  lie  both 
inside  and  outside  the  earth,  therefore,  the  general  expansion 
of  the  magnetic  potential  V  may  be  obtained  by  combining 
equations  (4),  (5)  and  (6b). 
'^   n 

V  =  a  \  N  Pg(e)[(cg(r/a)^  +  (l-cg)  (a/r)^-*"^}  Ag  cos  m0 
n=0  m=C 
+  {sn(r/a)''  +  (l-s^)  (a/r)''+^}Bg  sin  m^]  (10) 

v;here  eg  and  sg  are  numbers  ranging  between  0  and  1,  s/nd   rep- 
resent the  portions  of  the  harmonic  terms,  pg  cos  m0   or  Pg 
sin  m0,   in  V  which  at,  r  =  a,  are  due  to  matter  outside  the  • 
earth. 

The  observed  magnetic  elements  X,  Y,  and  Z  are  related  to 
V  by: 
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Y  =  -^v/r  sin  6^0) 

Z  =   av/^r 

Evaluating  these  expressions  at  the  earth's  surface,  r 

gives 

oo      n 

X    -  ^'  ^ 
•^r=a  -   / 


(11) 


=  a, 


d? 


m  r 


n_ 

J  de 


El 


An  cos  rx0   +  Bn  sin  m0 


(12) 


■r=:a 


n=:C  n2=0 


=  1/sin  e 


oo   n 

1  -c 1 


mPn 


An  sin  m0  -   b8  cos  m0 


(13) 


n=0  111=0 


oo 


n 


'r=a  - 


pS|'{ncg-(n+l)(l-c8)}Ag  cos  m0 


n=0  in=0 

+  {nsg-(n+l)(l-sS)}B8  sin  m0] 


(14) 


The  coefficients  Ag  and  Bn  of  equations  (12)  and  (I3)  are 
determined  by  a  series  expansion  of  Y  sin  6  and  X,  using  the 
measured  values  at  the  earth's  surface.   Details  of  this  nu- 
merical computation  are  given  by  Schffiidt  (1935) .   Since  An  and 
Bn  are  calculated  by  two  independent  means,  both  from  X  and  from 
Y,  disagreement  of  the  two  sets  of  values  indicates  a  non-poten- 
tial portion  of  the  earth's  field,  i.e.,  the  existence  of  a 
current  flow  transverse  to  the  earth's  surface.   The  reality  of 
such  current  systems  is  discussed  in  a  later  section  of  this 
report. 

Separation  of  the  geomagnetic  field  into  portions  of  ex- 
ternal and  internal  origin  proceeds  by  first  expanding  the 
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measured  values  of  Z  into  a  series. 


co      n 

"1  ^^^ — I 


Zr=a  =     }     /    '^niQ)i'^n  oos  m0   +^g  sin  m0)  (15) 

nr=0  ra=0 
Next,  by  comparinQ  equations  (14)  and  (15)  it  follows  that 
cKg  =  {ncg-(n+l)(l-cg)}Ag 
^g  =  {nsg-(n+l)(l-sg)}Bg  (16) 

IP         D. 

As  o^,  (i,  A,  B  are  knoim,  the  V3.1ues  of  en  and  sn  can  be  deter- 
mined . 

Gauss  (I839)  and  other  investigators  have  expressed  the 
magnetic  potential  in  terms  of  equations  (4),  (5),  and  (6b), 
where  the  coefficients  g^  ^'^'^-^   ^n  ^^®  termed  the  Gauss  coeffi- 
cients.  Such  coefficients  are  often  employed  to  compare  various 
analyses  performed  since  Gauss's  day.   Table  1,  due  to  Vestine 
(1947b),  compares  the  results  of  the  chief  analyses  made  of  the 
geomagnetic  field.   This  table  provides  a  rough  picture  of  the 
secular  variation  of  the  earth's  field  during  the  past  century, 
and  shows  that  the  major  part  of  the  field  is  due  to  a  dipole. 

Table  1.  The  first  eight  Gauss  coefficients  (external  plus 

internal  the  earth'  s  mar.netic  potential  V  expressed 
in  units  of  10  Y  . 


Source 

Date 

4 

si 

^i 

eE 

S2 

4 

A 

hi 

Gauss 

1335 

-3235 

-311 

■   625 

51 

292 

12 

-2 

157 

Erman  c 

1829 

-32C1 

-284 

601 

-8 

257 

-4 

-14 

146 

Peterson 

Adams 

1845 

-3219 

-278 

578 

0 

284 

-10 

4 

135 

Ada.ms 

18SC 

-3168 

-243 

603 

-49 

297 

-75 

61 

149 

Frits che 

1885 

-3164 

-241 

591 

-35 

286 

-75 

68 

142 

Schmidt 

1885 

-3168 

-222 

595 

-50 

278 

-71 

65 

149 

Dyson  6; 

1922 

-3C95 

-226 

592 

-59 

299 

-124 

144 

84 

Furner 

Afanasieva 

1945 

-3032 

-229 

590 

-125 

288 

-146 

150 

48 

Vestine 

1945 

-3C57 

-211 

581 

-127 

296 

-166 

164 

54 
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Transient  Variations  of  the  Earth's  Field 

A  discussion  of  the  transient  variations  of  the  geo- 
maGnetic  field  is  necessary  to  complete  a  description  of  the 
earth' s  magnetisn.   Such  variations  are  of  interest  since  they 
yield  valuable  information  about  the  behavior  of  the  sun,  the 
earth's  high  atmosphere,  and  outer  space.   Three  major  types  of 
variation  occur;  the  miaGnetic  disturbance  or  D  variation,  the 
solar  or  3  variation,  and  the  lunar  or  L  variation.  Each  mag- 
netic element  is  affected  in  a  characteristic  way  by  each  type 
of  variation. 

The  continuous  records  of  any  ms.gnetic  observatory  show 
that  on  some  days  all  magnetic  elements  exhibit  smooth  and 
re^jular  variations,  while  on  other  days  the  variations  are 
disturbed  and  quite  irregular.   Days  of  the  first  kind  are 
termed  quiet  days,  and  days  of  the  second  kind,  disturbed  days. 
For  each  observatorj''  a  figure  K,  ranging  geom-etrically  in  scale 
between  0  and  9,  is  assigned  to  indicate  the  value  of  magnetic 
variation  during  each  period  of  three  Greenv:ich  hours.   To  aid 
in  world-wide  analysis  of  magnetic  variation,  the  K   figure  is 
m.odlfied  to  a  iC^  figure,  v;hich  is  as  free  as  possible  of  the 
effects  of  local  m.agnetic  features,  and  is  dependent  upon  the 
latitude  of  the  station.   For  example,  Kp  =  9  would  indicate  a 
peak  magnetic  fluctuation  of  300  V  at  the  equator,  and  25C0Tat 
the  auroral  zones.  A  coarser  measurement  of  the  transient 
field  is  afforded  by  a  C  figure,  which  indicates  by  the  digits 
C,  1,  2  the  mean  value  of  the  magnetic  record  activity  for  each 
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GreenvTicli  day.  An  averac^e  of  the  C  figure  for  all  observ- 
atories provides  the  international  character  figure  for  that 
day.   Cn  the  basis  of  the  international  character  fi-gure  the 
five  quietest  days  (international  quiet  da^'-s),  and  the  five 
most  disturbed  days  (international  disturbed  da^'-s)  are  selected 
for  each  month. 

The  S  or  solar  daily  variation  is  a  function  of  local 
tinie,  season  of  the  year,  latitude,  and  sunspot  activity.   The 
variation  is  greater  and  more  rapid  during  daylight  hours  than 
during  darkness  hours.  A  spherical  harmonic  analysis  of  this 
variation  indicates  its  source  is  largely  external.  V.lien  de- 
termined by  averaging  hourly  values  of  variation  during  the 
five  quietest  days  of  the  month,  S  is  termed  Sq,  the  quiet-day 
solar  daily  variation.   In  considering  the  effect  of  S  upon 
each. magnetic  element,  LIcNish  (1939)  stated  that  the  ma>:imum 
varis.tion  of  X  is  at  the  equator  and  is  of  the  order  of  A0-6C^, 
the  maximum  variation  of  Y  is  at  50°  latitude  and  is  of  the 
order  of  35-50^,  a.nd  the  maximum  variation  of  Z  is  at  40° 
latitude  and  is  of  the  order  of  12-18^. 

The  lunar  or  L  variation  has  the  period  of  the  lunar  day, 
about  25  solar  hours.   Since  its  magnitude  is  only  about  1/15 
of  the  S  variation,  it  can  not  be  detected  directly  from  the 
magnetic  records,  but  must  be  calculated  by  statistical  methods, 
A  harmonic  analysis  shows  that  the  source  of  the  L  variation 
is  mostly  external. 

In  overall  behavior  the  D  or  magnetic  disturbance  varia- 
tion is  far  more  erratic  than  the  L  or  S  variations,  although 
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botli  a  27-d.ay  cycle  in  frequency  of  occurrence  and  a  general 
pattern  for  a  typical  disturbance  are  noted.   Kost  disturbances 
or  magnetic  storms  commence  at  all  points  over  the  earth  vrithin 
-|-  minute.   The  initial  phase  is  characterized  by  an  increase  in 
the  K  element  for  a  period  of  2  to  4  hours.   Then  H  falls 
greatly  below  the  undisturbed  value  held  before  the  Etorm. 
Lastly,  there  is  a  slov:  recovery  toward  the  normal  value  of  H, 
requiring  perhaps  several  days.   1-Iagnetic  declination  is  but 
slightly  affected  during  the  storm  period;  the  changes  in  ver- 
tical intensity  Z  are  much  smaller  than  the  changes  in  H,  a.nd 
have  an  opposite  sign.  In  several  severe  magnetic  storms  an 
H  variation  up  to  2000^  has  been  noted. 

G-eographically,  the  magnitude  of  the  D  disturbance  in- 

o 
creases  from  low  latitudes  to  a  maximum  at  65  magnetic  lati- 
tude, the  auroral  zone,  followed  by  a  slight  decrease  as  the 
poles  are  approached.  At  liigli  latitudes  few  days  are  com- 
pletely free  of  some  magnetic  disturbance.   'The  D  variation  is 
connected  with  solar  activity?-,  but  a  time  lag  of  one  to  three 
days  occurs  betv;een  peal-:  sun-spot  activity  and  the  onset  of  a 
period  of  terrestrial  magnetic  storms. 

MODELS  CF  TEE   EARTH'S  FIELD 

2<Iain  Dipole  Field 

In  order  to  first  develop  a  model  of  the  earth's  field, 
a  comparison  is  miade  between  the  various  components  of  the. 
field  due  to  internal  and  external  sources.   Table  2  lists  the 
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internal,  external,  and  conbined  Gauss  coefficients  of  tvo 
analyses  made  from  r.aGnetic  survey  data  for  the  year  19^5 
(Vestine,  196Ca;  Fa.nselau  and  Eautzleban,  1958) .  As  the  co- 
efficients relating  to  an  external  field  source  arc  very  snail, 
it  is  uncertain  that  the  earth's  lon^-term  field  has  an^""  ex- 
ternal sources.   Only  the  transient  cagnetic  variations  D,  S, 
and  L  appear  to  hr^ve  external  sources. 

Table  2.   Values  of  Cn  =  f(Sn)^  +  (hn)^}"^'  for  field  of  origin 
external  and  internal  to  the  earth,  units  of  10 ^'^ 

Source 

Fanselau-Kautzeben 
n         0 

u  It 

Vestine 
n 

Line  integrals  of  H  a.bout  closed  paths  on  the  ee.rth'  s 
surface  have  been  shoi-m  tc  be  non-zero  with  randonly  distrib- 
uted values.  All  the  Eieasurements  of  earth-air  currents  normal 
tc  the  earth's  surface  are  far  toe  small  to  account  for  the 
non- vanishing  of  these  integrals.   It  seems  lihely,  therefore, 
that  if  the  magnetic  survey  data  possessed  greater  accuracy  and 
more  detail,  these  line  integrals  would  have  zero  values.   If  ' 
the  results  of  this  line  integration  do  have  physical  reality, 
the  portion  of  the  field  attributable  to  a  vector  potential 
woLild  be  less  than  1/^  of  the  entire  field.  No  serious  error  is 
introduced  by  assuming,  as  in  the  spherical  harmonic  analysis, 
that  the  earth's  surface  field  is  derivable  entirelj''  from  a 
scalar  potential. 

The  first  order  terms  of  the  me.gnetic  potential  may  be 
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616 

340 
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Inner 
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Inner 
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122 
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I'/ritten  for  r  =  a  as 

V/a  =r  Qi   cos  e  -5-  si   cos  0  +  li^  sin  0^   sin  8.         (1?) 
Now  Hq,  ©O'  00   ^^®  defined  as 

cos  00  =  g5/Ho 

tan  00  =  lii/sj 
The  a.n^le  0  betv.-een  the  direction  (0,0)  and  the  direction 
(©O,0o)  ^°   expressed  "by 

cos  0  =  cos  e  cos  Qq   +  sin  0  sin  Qq   cos(0-j^o)        (^5) 

The  iiia£pietic  potential  nov?  "becomes 

V/a  =  Ko  cos  0  (20) 

This  potential  at  the  earth's  surface  cculd  "be  the  result 
of  a  magnetic  dipole  moment  Hoa^  s.t   the  center  of  the  earth, 
with  direction  (-0o>  -0o)  •   ^^^-®  s.xis  of  this  dipole  intersects 
the  earth's  surface  at  points  termed  the  geomasnetic  poles.   At 
these  poles  the  vertical  magnetic  intensity  due  to  the  dipole 
will  be  ±2Eqj   and  at  90°  from  these  poles  the  horizontal  inten- 
sity will  be  Hq.   If  the  earth  were  uniform.lj'-  magnetized 
throughout  v;-ith  a  magnetization  of  3Ho/4Tr  cgs  poles  per  cm.  , 
a  centered  dipole  field  would  also  result. 

A  third  miodel  is  a  surface  current  S3^stem  flowing  east  to 
west  and  having  a  distribution  (Chapman  and  Bartels,  19A0), 

I  =  30/4 TT  aKo  cos  0  amperes  (21) 

or  a  current  density 

J  =  -?l/a.dQ   =  0.8  sin  0  amperes/cm.  (22) 

An  infinitude  of  other  models  can  also  TDroduce  the  first  order 
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field  observed  at  the  earth's  surfe.ce,  out  the  three  preceding 
ncc.els  are  relatively  cinple  to  visualize. 

iLn  analj'-sis  of  the  Geomagnetic  field  by  Finch  and  Leaton 
(1957)  J  using  data  for  the  j-ear  1955,  Qlves   tlie  value  of  Hq  as 
0.3120  Qcxisr^,    and  the  ma^-.netic  moment  of  the  earth's  centered 
dipole  as  8.O6A  X  10^  cgs  pole  cm.  Hq  has  been  decreasing  at 
the  rate  of  15  ^ a  year  during  the  past  century,  although  recent 
surveys  indicate  that  this  rate  of  decrease  may  noi;  be  ap- 
proaching zero.   The  geomagnetic  poles  are  located  presently 
at  79°N,  290°E,  and  79°S,  llO^S.   These  poles  should  not  be 
con-fused  i;ith  the  ms,gnetic  dip  poles,  which  are  points  of  zero 
vertical  intensity  Z,  and  are  located  asjrmmetrically  at  about 
73°N,  2620E,  and  68^S,  145°E.   It  is  not  certain,  judging  from 
analj'-ses  telien  over  the  past  century,  that  the  geomagnetic 
poles  are  shifting.  However,  the  dip  poles  var^^  rapidly  in 
position  during  an  otherwise  sliglit  change  in  the  geomagnetic 
field  pattern;  their  position  has  little  significance  in  an 
analysis  of  the  main  field. 

In  the  study  of  such  phenomena  as  magnetic  storms,  cosmic 
rays,  or  auroras  it  is  often  convenient  to  adopt  a  geom.agnetic 
co-ordinate  system  tahen  with  reference  tc  the  a::is  of  the   cen- 
tered dipole.  I-IcNish  (193^)  ^''-as  developed  a  set  of  nomographs 
that  give  the  geomagnetic  co-ordinates  for  any  point. 

If  the  position  of  the  main  field  dipole  is  cliifted  wliile 
keeping  its  vector  moment  constant.  Chapman  and  Bartels  (19^0) 
have  shov-TL  that  tlie  terms  in  V  involving  P?  ^nd  P2  may  be  made 
to  vanish.  Such  a  model  is  that  of  an  eccentric  dipole;  its 
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present  position  is  3CC  km.  from  the  earth's  center  in  the  di- 
rection 15°K,  168^2.   This  eccentric  dipole  fits  the  magnetic 
data  of  the  Southern  hemisphere  somewhat  better  than  does  the 
centered  dipole;  there  is  little  difference  betv/een  the  tv/o 
models  over  the  Northern  hemisphere. 

Residual  Field 

If  the  potential  of  the  centered  dipole  is  subtracted  from 
the  total  potential  of  the  geomagnetic  field,  the  remaining 
potential,  containing  spherical  harmonic  terns  of  degree  two 
and  higher,  is  that  of  the  residual  field.   An  iscdynamic  chart 
of  any  residual  field  element  has  irregularly  distributed  re- 
gions of  high  and  low  field  intensity  v;hich  range  in  diameter 
from  25*^  to  10C°.  Maps  of  the  secular  variation  of  the  earth's 
field  also  shov;  broad  regions  of  increasing  and  decreasing  field 
v;ith  dimensions  quite  similar  to  the  distribution  of  the  re- 
sidual field.   It  can  also  be  demionstrated  that  the  features 
of  both  the  non-dipole  field  and  the  secular  variation  have 
drifted  v/estward  during  the  past  half  century.   Thus  there  is 
every  reason  to  believe  that  the  tv70  phenomena  are  correlated 
such  that  the  secular  variation  is  the  time  derivative  of  the 
residual  field. 

Vestine  (lS47a,b)  published  a  detailed  spherical  analysis 
of  the  earth's  field  for  the  period  1905  to  1945.  Using  the 
data  of  Vestine,  Bullard  et  al.  (195C)  computed  the  angular 
shift  of  the  field  for  1905  that  would  best  fit  this  field  to 
the  field  of  1945.   As  a  result,  the  drift  of  the  non-dipole 
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field  was  riven  as  CIS  wectvarci  -per  jeer,    and  the   drift  of 
the  secular  variation  as  0.32^  westward  per  yee.r.      Both  drifts 
were  essentially  independent  of  latitude.   It  is  not  locio\rn 
wlietlier  these  drifts  can  he  extrapolated  into  the  past,  as  the 
older  r.iaQ;netic  survey  data  was  not  sufficiently  complete  to 
develop  an  accurate  nap  of  the  residual  field.  However,  El- 
sasser  (1950a)  pointed  out  that  the  intersection  of  the  west- 
erly line  of  zero  magnetic  declination  with  the  equator  he.d 
shifted  90°  westward  in  AOC  j'-ears'.   Since  other  features  of  the 
older  survey  data  also  support  the  possihilit3^  of  a  westwe.rd 
drift  of  the  non-dipole  field,  this  drift  is  now  s.ccepted  as  a 
physical  reality.  Therefore,  any  theory  about  the  crici^-  of 
the  r;eonetric  field  must  adequately  explain  thj.s  drift. 

Seismic  data  has  indicated  that  the  earth  has  a  core  of 
dense  fluid,  probably  molten  iron  and  nichel,  with  a  radius  of' 
about  O.A-5   earth  radii.   Thus  most  models  of  the  earth's  re- 
sidual field  have  placed  the  sources  at  the  surface  of  the 
molten  core,  since  any  time  varying  field  in  the  interior  of 
the  core  would  be  screened  out  'oir   the  high  conductivity  of  the 
core.   !-IcIIish  (1940)  represented  each  nest  of  vertical  inten- 
sity isodynar.ic  linos  b3'-  a  radially  directed  dipole  placed  at 
r  =  a/2.  Vestine  (1947b)  assumed  current  fimcticns  of  the  form. 

In  =  10  2n-^l  V^  (a/r)^"^^  (23) 

A- TV        n 

where  Vp_  is  the  potential  at  the  earth'  s  surface,  and  In  is  the 

current  distribution  in  a  thin  current  sheet  at  the  radius  r. 

The  total  current  function  that  would  generate  the  residual 
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field  vras  t-Iien  calculated  and  r:raplied  for  current  cheet  depths 
of  0,  lOCC,  20CC,  and  30OO  Ian.   Lowes  and  Runcorn  (1951)  used  a 
Graphical  analysis  on  the  X  and  Y  ccnponents  of  the  field  to 
obtain  a  dipole  distribution  fitting  the  observed  secular  vari- 
ation field  to  i'.-ithj.n  1^   accurac3''.   Ivjelve  vertical  dipoles. 
placed  just  beneath  the  surface  of  the  earth' s  core  were  neces- 
sary; each  dipole  represented  a  horizontal  circuit  of  current 
flow  at  the  surface  of  the  core. 

External  Field 

A  complete  model  of  the  rz^eoiLO.ryietlc,   field  nust  include  not 
only  the  earth's  internal  field  for  v/hich  models  have  now  'DQon 
given,  but  also  tlie  transient  maj^netic  variations.   Only  models 
for  the  So  variation  and  the  D  variation  are  considered  here, 
as  these  two  variations  provide  the  major  portion  of  the  tran- 
sient fluctuations.   Spherical  harmonic  analysis  indicates  that 
the  source  of  both  tlie  Sq  sjid  D  variations  is  a.bout  75?^  ex- 
ternal and  25/^  internal  to  the  earth.   The  internal  portion  of 
both  variations  is  due  to  currents  induced  by  the  external 
field  in  the  mantle  of  the  earth.   Just  as  in  the  case  of  the 
earth's  main  field,  there  are  innum.erable  models  that  would 
reproduce  the  observed  transient  variations.   This  discussion 
will  consider  only  those  models  that  a-re  consistent  with  the 
present  Imowledse  of  the  transient  variations. 

The  quiet  ds-y  solar  variation,  Sq,  is  usually  considered 
to  have  its  source  in  the  ionosphere  at  a  height  of  70  to  400 
kn.   From  observations  of  S^  Chapman  and  Bartels  (19^0)  have 
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developed  an  earoli  wide  ovcrliead  current  system  that  v/culd 
produce  tills  fluctuation  at  tlie  earth's  surface.   This  current 
s^'ston,  V'hich  neves  v:esti/ard  vrith  the  sun  as  the  earth  rotates, 
consists  of  four  circuits  cr  nests  of  flo\r   lines;  t'.:o  circuits 
are  in  the  horthern  hemisphere,  ?aid  t'-ro  are  in  tlie  Southern 
hemisphere.   In  either  hemisphere  one  circuit  covers  the  day- 
light half  cf  the  earth;  the  second  circuit,  carrji^in;;;:  less  cur- 
rent thaj.1  the  first,  covers  the  darh  half  of  the  rlohe.   The 
overall  structure  of  the  current  s^^stem  varies  considerably 
■'.v'ith  the  seasons;  the  circuits  shift  scuthvjard  during  summer  in 
the  Northern  hemisphere  and  northward  during  summer  in  the 
Southern  hemisphere . 

The  effect  cf  the  mia^gnetic  disturbance  variation  D  can  be 
Gim.ulated  by  two  superimposed  current  systems.   Cne  portion  of 
D,  term.ed  the  storm- time  variation  Bgrpj  is  a  function  of  the 
time  elapsed  since  the  beginning  of  the  disturbance,  and  ca/n  be 
represented  by  an  overhead  current  system  that  is  symmietric 
about  the  earth's  axis.   Cliapman  and  Bartels  (I9^f0)  plotted 
such  an  idealized  current  system  for  an  arbitrary  height  in  the 
atm.osphere.   Their  current  system  for  D^^   has  its  greatest  cur- 
rent density''  at  the  equator,  with  decreasing  current  densities 
e.Q   the  poles  are  approached,  illfv^n  (1955),  2,nd  also  Appel 
et  a,l.  (1962)  postulated,  as  a  m.ore  probable  source  of  the  Ds? 
variation,  an  equatorial  ring  current  girdling  the  earth  at  a 
distance  cf  several  earth  radii.  I-Iea.Gurements  made  hj   both 
Russian  and  American  satellites  tend  to  verify  that  several 
such  ring  ciirrents  exist,  one  in  the  outer  Van  Allen  belt, 
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another  ■'.•-ell  be^^ond  the  Van  Allen  "belt  at  about  6  to  10  earth 
radii  (Ahasif\i  and  Chapman,  1961;  Smith  et  al.,  I960). 

The  second  portion  of  tlie  D  variation  is  Sj),  a  daily 
fluctuation  that  occurs  only  during  periods  of  na^notic  storm 
activit^^.   The  current  systen  responsible  for  Sj)  seenis  to  be 
located  in  the  ionosphere;  its  position  is  asymmetric  \ilth   re- 
spect to  the  ea.rth' s  axis,  G.nd  follov:s  the  sun  as  the  earth  ro- 
tates. At  the  polar  recions  the  Sj)  current  system  has  its 
higiiest  current  density.  Over  the  earth  its  currents  floii   east 
and  west  along  the  equator,  and  north  and  south  alcn{^  the  noon 
meridian.  An  idee^lized  current  system  for  D,  combining  the  tuo 
current  systems  for  Dst  Q-^^  Sp,  has  been  developod  ''oy   Cliapman 
and  Bartels  (19^0). 

TI-IECRIES  OF  THE  EARTH'S  INTERl^AL  FIELD 

Es,rly  Theories 

In  spite  of  the  ma-n^-  theories  devised  in  the  past,  the- 
source  of  the  earth's  internal  field  is  still  a  highly  specu- 
lative topic.  Any  satisfactory  theory  must  explain  the  field's 
predominantly  dipolar  chare-cter,  the  orientatic^n  2.nd  rate  of 
cliange  of  this  dipole  field,  the  rapidly  varying  behavior  of 
the  non-dipole  field,  and  the  earth's  past  ma,gnetic  historj''  ob- 
tained from  fossil  magnetism.  Since  the  inaccessibility  of  the 
earth' s  interior  to  direct  observation  prohibits  a  conclusive 
test  of  cAij   theory,  the  various  theories  are  compared  in  order 
to  select  the  hyi'iothesis  r-ost  closely  predicting  the  presently  ' 


Imouii  GGomaryietic  plienomena. 

Possibly  tliG  first  explanation  of  the  origin  of  the  in- 
ternal field  uas  the  assunption  of  sono  arbitrary  distribution 
of  forromaEnetlc  inateric.l  irithin  the  earth.   Kovjever,  the  In- 
crea.se  of  the  eartli'  s  ter.iperature  vrith  depth  seems  to  indicate 
that  the  Curie  point  of  nost  ferrona^rnetic  naterials  would  be 
reached  at  a  depth  of  25  hn.   Since  a  Curie  point  increase  with 
increa-sing  pressure  ha,s  never  been  verified  experimentally,  any 
distribution  of  nacnetic  ns.terial  nust  be  restricted  to  the 
cnistal  region .  A  calculation  of  the  necessary  a-verae©  mag- 
netization in  the  crust  gave  the  valve  of  5  to  8  cgs  poles  per 
cm.  ,  a  value  more  than  twice  the  naxiinum  remanent  magnetiza- 
tion of  iron  or  nichel,  and  over  100  times  the  permanent  mag- 
netism of  igneous  roclis.  Besides  requiring  these  impossible 
values  of  magnetization,  the  ferromo^gnetic  theory  could  not  ex- 
plain the  rapid  variations  of  the  non-dipole  field. 

ilnother  earlj^  theory  postulated  that  the  ms.gnetic  field 
was  produced  by  the  rotation,  with  the  earth,  of  a  system  of 
electric  charges.   This  system,  could  be  either  a  uniform  sur- 
face cha.rge  S,  or  a  volume  charge  Q  uniformly  distributed 
tliroughout  the  earth.  However,  the  earth's  conductivity  seoms 
great  enough  to  destroy  opaite  rapidly  any  such  charge  distri- 
bution. 

Blackett  (1947),  noting  the  existence  of  axial  m-agnetic 
fields  in  various  rotating  cosmic  bodies  other  tlie-n  the  earth, 
proposed  a  universal  relationship  between  angtilar  momicntum  and 
magnetic  moment  in  large  bodies;  any  rotating  body  would  thus 


27 


be  spontaneously'  inagnstized.   For  the  case  of  the  earth  such 
a  funda^mental  hypothesis  cculd  he  directly  tested;  the  source 
of  the  field  v;culd  he  distributed  uniformly''  throughout  the 
whole  earth,  and  the  vertical  component  of  field  intensity 
would  not  obey  an  inverse  cube  law  viithin  the  earth's  crust. 
Measurements  of  vertical  intensity  carried  out  in  coal  mines 
by  Runcorn  et  al.  (1951)  c-i^  not  support  this  theory,  and 
Blackett  (1951)  soon  announced  the  results  of  a  "negative" 
laboratory  experiment  upon  his  theory. 

As  the  time  constant  of  a  system  of  freely  deca3''ing  cur- 
rents wiohin  the  earth  is  estims-ted  to  be  about  2  million  years, 
such  a  system  might  seen  to  offer  an  answer  to  the  problem  of 
the  geomagnetic  field.   ITith  the  assumption  that  the  current 
system,  was  started  100  m.illion  years  ago,  a  slicrt  interval  of 
timie  compared  v:ith  the  probable  age  of  the  earth,  the  required 
initial  current  would  have  been  about  10  '  tiroes  the  presently 
required  value  of  lOCO  million  am.peres.   No  process  is  knov/n 
that  could  set  up  such  enormous  currents  in  the  earth. 

Nevertheless,  internal  current  systems  s.re  a  llliely   source 
of  the  earth's  field,  but  some  m.echanism  is  required  to  m.ain- 
tain  these  systemis.   One  proposed  m.echanism  is  the  thermoelec- 
tric effect;  a  second  mechanism  is  electromagnetic  induction. 
Both  mechanisms  require  the  existence  of  ccnvective  tlovr   pat- 
terns in  the  liquid  core  of  the  earth.   The  existence  of  such 
fluid  motions  is  reasonably  certa^in,  for  only  a  small  amount  of 
heat  is  necessary  to  initiate  convection;  this  heat  could  de- 
rive from  radioactive  heating  or  slovr  chemical  se-oaratlons  and 
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crystallizations.   The  rapidly  varying  nature  cf  the  ncn-dipole. 
field  also  points  cut  the  physicG,l  realit^^  of  cor^vectlve  flov; 
patterns  in  the  earth's  core;  the  fluid  motions  interact  v/ith 
a  current  system  to  produce  the  irregular  distribution  of  the 
secular  variation  and  residual  fields.   There  is  more  than  a 
superficial  resemblance  betv.-een  v/eather  maps  a.nd  magnetic 
charts;  both  display  the  phenomena  cf  large  scale  fluid  motion. 

Elsasser  (1939)  su^-^ested  that  thermoelectric  potentials 
were  the  source  of  an  internal  current  system  producing  the 
Seomiarnetic  field.   This  concept,  extended  by  Huncorn  (195^), 
has  som.e  irxiierent  difficulties  that  mahe  it  less  satisfactory 
than  the  induction  or  dynamo  theory.  A  thermoelectric  e.m.f.  is 
developed  at  the  surface  of  the  core,  the  core  and  mantle  form- 
ing the  tv/o  materials  of  the  thermo Junction.   Convection  cur- 
rents in  the  core  must  miainta,in  temperature  differences  between 
various  points  on  the  boundary  of  the  core.   The  resultant  cur- 
rent system  generates  a  magnetic  dipole  moment  only  through  the 
action  of  the  Coriolis  force  on  the  convecting  fluid  streams. 
It  is  shovrn  in  the  discussion  on  magneto-hydrod^niamdcs  that  a 
magnetic  line  of  force  contained  in  a  conducting  fluid  acts  as 
if  it  were  attached  to  one  set  of  fluid  particles,   vrnen  the 
fluid  moves,  the  line  of  force  is  deformed  in  following  the  mo- 
tion cf  the  fluid  particles.   Thus  the  Ccriolis  force  in  de- 
flecting the  fluid  flow  distorts  the  em.bedded  m.agnetic  field' 
and  produces  a  net  magnetic  dipole  m-oment. 

One  difficulty  cf  the  thermoelectric  theory  is  that  it 
must  include  one  of  the  follov.'ing  requirem.ents:  (1)  large  • 
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tenperature  differences  in  tlie  core,  (2)  lar^e  tlicrno electric 
coefficients,  or  (3)  extreme  georn.etry  of '  ccnvective  flow. 
Also,  the  sic'.n  of  the  thernoelectric  coefficients  is  ur-laipvrn, 
and  no  explanation  of  the  probable  past  reversals  of  the  ea.rth'  s 
marpietic  polaritj'  is  ^iven.   There  is  sone  evidence  (Vestine, 
195^)  th-C-t,  dv.e  to  hJ.rrh  pressure,  the  nantle  material  is  a  seni- 
conductcr  and  possesses  a  reasonabl^r  hi^h  Hall  coefficient. 
Vestine  then  5U26®sted  that  the  thermoelectric  and  Hall  effects 
interact  in  a  regenerative  process  responsible  for  at  least  a 
part  of  the  earth's  internal  field. 

Basic  !-a£7ieto-hydrod3rnamics 

Before  proceodin^  to  a  detailed  discussion  of  the  d;^^nano 
thqory,  it  is  desirable  to  consider  the  basic  eouations  of 
ma^^^e  to -hydrodynamics.  Any   internal  dj'n.amo  action  must  satisfy 
these  equations. 

If  the  displacement  current,  all  purely  electrostatic 
effects,  and  relativistic  effects  of  order  higher  than  v/c  are 
Ignored,  lls.xwell' s  equations  in  the  rationalized  ml^s  system. 
are : 

curl  E  =  -BB/at  (24) 

curl  3  =^J  _  (25) 

div  B  =  0  (26) 

i;here  E,  B,  and  J  are  the  vectors  of  the  electric  field,  Tia^- 
netic  flux  density,  and  current  densitj'-,  respectively.   The 
ns.rp.Qtlo   permeability  is^  . 

For  an  electrically  conducting  fluid  having  relative 
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motion  v:ith  respect  to  the  rr.agnetic  field,  the   current  density 
is  • 

J  =  cj^(E  +  V  X  B)  (27) 

with  d'  as  the  electrical  conductivity,  and  v  as  the  velocity 
vector;  dv   X  3  boins  the  motional  induction  term. 

Taking  the  curl  of  equation  (25),  and  usino;  equations  (24) 
and  (27)  to  eliminate  E  gives 

curl  curl  B  =  Xf^  -^B/^t  +  curl  (v  X  B)  (28) 

Nov;, 

curl  curl  B  =  crad  div  B  -  ^B  (29) 

where  X^  3  is  the  vector  Laplacian  of  B.  Combining;  equations 
(26),  (28)  and  (29)  yields 

^3/at  =   curl  (v  X  B)  +  Vj^  2^1  (30) 

The  quanity  l/x/a^  has  been  ^■.Tr'itten  as  v^  and  is  called  the  mag- 
netic viscosity. 

The  Hapier- Stokes  equation  of  hydrodynamics  for  a 
moving,  electrically  conducting  fluid  is 

.  av/^t  +  (v  .V)v  =  -1/pVp  -VU  +  vS^v  +  1/p  J  X  B    (31) 
with  p  the  fluid  density,  U  the  gravitational  potential,  v  the 
kinematic  viscosity,  and  p  the  pressure.   The  term,  l/p  J  X  E, 
is  the  pondermotive  force  v/hich  the  electromagnetic  field  ex-  . 
erts  on  the  fluid,  and  represents  the  coupling  between  the 
field  and  the  fluid.  A  second  necessary  equation  of  hydro - 
dyna-mics  is  the  continuity  equation  which,  for  an  incompressible 
fluid,  is  ■  •        • 

div  V  r:  0  (32) 

Using  equation  (25)  in  equation  (3I)  results  in 
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^v/^t  -^-  (v  .V)y  =  -1/.^  V  p  -VU  +  V  Vv 

+  1/p/^  curl  B  X  3  (33) 

Equations  (3C,  (32),  and  (33)  contain  onl^r  the  vectors 
V  and  3,  and  are  the  basic  equations  of  field  motion.   The 
equations  include  non-linear  terms  of  three  kinds:  the  term 
curl  (v  X  E),  the  term  1/qju    curl  3X3,  and  the  inertia  term 
(v  rV)v.   ihese  non-linearities  make  a  mather.atical  solution 
almost  impossible.   Houever,  bj.  a  dimensional  analysis  Slsasser 
(I95^a)  obtained  scr;e  important  c.ualitative  results. 

This  analysis  begins  by  combining  equations  (25)  and  (27) 
into 

E  +  V  X  B  =  \/2i  curl  3    •       ,  (34) 

If  L,  T,  and  V  represent  the  orders  of  magnitude  of  length, 
time,  and  velocity'',  respectively,  then,  03^  equation  (24),  S  has 
the  order  of  magnitude  VE,  as  has  the  v  X  3  term.,  wliile  the 
term.  Vr-  curl  3  is  of  the  order  u-^-T'^B, 

In  hydrodynamics  the  nondimensional  Reynolds  number  is 
define-d  as  Pw  =  IM/v,      Thus  the  viscous  term  v^v  in  equation 
(33)  is  o-  the  order  1/R  compared  to  the  inertia  term  (v  •\7)v. 
A  m.agnetic  Reynolds  number  may  be  defined  as 

R^  =  LV/^p^  (35)' 

The  terms  magnetic  viscosity  and  magnetic  Re3rnolds  number  can 
be  justified  by  noting  that  the  dissipative  term.  \/_^^3  is  of 
order  l/Rn  compared  to  the  other  terms  in  equation  (3C)  a  re- 
lation analogous  to  the  purely  hydrodynamic  case. 

The  right-hand  side  of  equation  (34)  is  of  the  order  l/R^i 
•v-:hen  compared  to  the  left-hand  side.  Tvlien  Rja  takes  on  large 
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values  due  to  L  becoming  lar^e,  as  V,  B,  and.  v^   renain 
constant,  equation  (3^)  ^ives 

E  =  -V  X  B  (36) 

Therefore,  the  actual  electric  currents  corresponding^;  to  a 
given  nagnetic  field  can  "oGccrie  very  sr.all  in  systems  of  cosmic 
dimensions. 

A  sir-ple  physical  interpretation  of  R^   is  gained  by  consid- 
ering an  electric  conductor  in  v/hich  a  system  of  currents  floivs. 
If  electromotive  forces  are  absent,  these  currents  will  decay 
exponentiall3'^.   The  decay  time  Z   can  be  estimated  from  equation 
(30)  hy   setting  v  =  C,  and  is  of  the  order 

T  =  L-/v/^  (37) 

so  that 

%  =  -CV/L  =  X/T  (36) 

where  T  is  the  time  for  a  fluid  particle  to  tr>avel  a  distance 
equivalent  to  the  extension  L  of  the  conductor.   If,  as  in  the 
case  of  the  earth's  iPiterior,  Rj^  has  a  l8.rge  value,  then  the 
decay  timetsiust  be  very  large  compared  to  a  period  of  Fxechan- 
ical  motion.   Fluid  mictions  may  be  considerable  during  a  time 
in  V7hich  the  decay  of  the  electromiagnetie  field  is  small. 

One  more  general  result  of  m.agneto-hydrodynamics  is  found 
by  com.bining  equations  (24)  and  (30),  and  v:riting  them,  in 
integral  form  as 

.^3/^t  .•  dA  =  j  V  X  B  .•  dl  -  Vjjj  /curl  3  ?  dl       (39) 
Next  let  a  surface  S  and  its  contour  miove  bodily  with  the 
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fluid,  sweeping  out  a  disk-shaped  volume  during  time  dt.  The  ■ 
flux  passing  through  the  strip  that  forms  the  edge  of  the  disk 
will  be,  by  Fig.  3, 


d$  =  /  /  (dr  X  dl)  .  B  =  - 


V  X  B  •  dl 


dt  (40) 


Thus  the  total  time  derivative  of  flux  crossing  S  as  it  moves 
with  the  fluid  will  be  -       • 

d/dty  1  B  •  dA  =  ^^^^J  J  B.dA-/vXB-di        (4l) 

The  first  term  on  the  right  is  the  flux  changing  with  time  for 
a  surface  fixed  in  space;  the  second  integral  on  the  right 
represents  the  change  in  flux  for  a  fixed  B  as  the  surface  is 
displaced.   Substituting  equation  (41)  into  equation  (39) 
provides 

d/dt  /  /  B  •  dl  =  -v/in  /  curl  B  •  dl  (42) 

For  ideal  conductors  1/^=0  and 

d/dt  /  /  B  .  dl  =  0  (43) 

Equations  (42)  and  (43)  hold  equal Ijr  v;ell  for  both  in- 
com.pressible  and  com.pressible  fluids,  and  indicate  that  the 
magnetic  lines  of  force  in  a  conducting  fluid  are  carried  bodily 
along  with  the  fluid  motion.  .  For  a  fluid  of  finite  conductivity 
the  right-hand  side  of  equation  (42)  represents  a  "slippage"  of 
the  lines  of  force  when  follov/ing  the  fluid  motion. 


3^ 


s(t) 


S(t  +  dt) 


FIG.     3 . 
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The  Dynamo  Theory 

The  djrnaxo  theory  has  the  widest  acceptance  of  an^^  theory 
of  the  earth's  interno.l  field  since  it  presupposes  the  least  in 
the  i:ay  of  natural  phencLiena  and  can  easily  predict  much  of  the 
behavior  of  the  field.   The  equations  of  magneto-hydrodynamics 
suGS^st  that  the  magnetic  field  is  produced  "by  the  mechanism  of 
electroriagnetic  induction;  magnetic  energy  is  drai-ni  from  the 
kinetic  energy  of  the  fluid  in  the  core.   Thus,  B  and  v  fcrmi  a 
single  dynamic  system;  energy  is  not  separately  conserved  for 
each  of  these  fields. 

Since  the  ba,sic  equations  involving  B  and  v  are  non-lines^r, 
it  is  not  possible  to  obtain  anal^/tically  a  steady-state  d^iiamo 
solution  satisfying  the  necessary  boundary  conditions.   All  ap- 
proaches taken  so  far  have  separated  the  hydrodynamic  and  the 
electromagnetic  problems  and  attempted  to  solve  only  the  latter 
problem.   Some  velocity  distribution  ie  assumed,  and  inter- 
actions betv:een  the  magnetic  field  and  the  fluid  motion,  satis- 
f3''ing  Kaxvjell's  equations  and  capable  of  maintaining  the  field, 
are  scuglit. 

An  im.portant  characteristic  of  any  mechardcal  d3''namo  is, 
its  low  order  of  symmetry;  the  armature  has  the  S3amrietry  of  a 
clock  face  in  which  the  two  directions  of  rotation  are  not 
equivalent.   Depending  upon  the  rotational  direction,  such  a 
dynamo  process  will  be  either  degenerative  or  regenerative.  A 
simple  body  like  a  sphere  has  a  high  order  of  structural  sym- 
metry.  Thus  to  act  as  a  djiiamo  the  earth's  core  m.ust  possess 
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some  as;^ininetr3^  in  fluid  motion,  for  Cowling  (193A)  hc.s  shovm 
that  hlrrhly  s^^rjnetrical  fluid  notions  in  a  sphere  carjnot  pro- 
duce e.ny  dynano  action.   The  necessary  asyiisHetry  in  fluid  flovr 
is  most  likely  provided  by  the  Coriclis  force.   As  in  the 
mechanical  djmamo,  the  resultant  induction  reaction  between 
the  fluid  and  the  field  must  have  the  proper  si^n  in  order  that 
the  field  be  maintained. 

For  a  sphere  the  vector  fields  of  E,  B,  J,  and  v  can  be 
derived  from  the  scalar  functions: 
y  gs=  R(r)  P^  m^cos  e)  sin  :n0 

The  fields  derived  from  these  functions  are  of  two  hincs:  a 
toroidSvl  field  T  =  curl  (rV)  =  grad  V  X  r,  and  a  poloidal  field 
S  =  curl  T. 

To  produce  the  poloidal  field  S-j  (corresponding  to  the 
main  dipole  field),  Bullard  and  G-ellman  (195Aa)  postulated  thr^t 
a  much  more  powerful  toroidal  magnetic  field  must  exist  in  the 
earth's  core.   This  toroidal  field  would  not  appear  in  any 
measurements  made  at  the  earth's  surface.  Bullard  also  shewed 
that  no  single  type  of  velocity  field  could  provide  d^mamo 

action.   Of  the  combinations  of  two  velocity  fields  possibly 

2c 

allowing  d3''namo  action,  the  process  requiring  !]_  and  S2   fluid 

2c 

m-otions  was  the  simplest.   The  details  of  this  T-j_S2  inter- 
action are  given  in  Fig.  4  and  Table  3.  After  extended  numer- 
ical  computation,  it  was  concluded  that  the  T^So  interaction 
could  produce  d^mamio  action.   This  particular  solution  is  ^y   no 
means  unique,  and  somiO  doubt  ha,s  been  cast  upon  the  convergence 
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Coupling  by    Tj 
motion 

Coupling   by  S^^ 
motion 


FIG.  4.     RELATIONS      BETWEEN     INDUCED    AND    INDUCING     FIELDS 
FOR     A     1|S|^     DYNAMO     PROCESS. 


Ta'Dle  3.  .   Relations  betv;een  induced  and  inducing 
fields   for  Ti    and  S§*^   fluid  motions. 
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of  the  solution.   Nevertheless,  Bulle/rcL's  results  ir^dicate  that 
the  internal  dj'-nano '  theory  is  quite  feasible  and  can  provide  a 
self-consistent  account  of  the  earth' s  magnetic  field. 

In  order  to  shov:  that  a  marr.netic  field  could  be  defcrned 
in  such  a  way  so  as  to  rroduce  a  self-sustaining  dynamo,  Parker 
(1955)  and  Elsasser  (1955)  considered  two  effects  of  the  earth's 
rotation  upon  fluid  notions  in  the  cere.   First,  a  simple  Tj 
toroidal  magnetic  field  -.-Tith  lines  of  force  circling  about  the 
axis  of  the  earth,  parallel  to  lines  of  latitude,  is  assumed. 
Ihis  toroidal  field  results  from,  the  deformation  of  the  poloidal 
S-]_  (dipole)  field  by  the  non-uniform  rotation  of  fluid  particles 
in  the  earth's  cere  (Fig.  5).   The  toroidal  field  cs.n  be  gener- 
ated directly  from  the  poloidal  field,  but  not  vice  versa. 

The  production  of  a  poloidal  field  from  the  toroidal  field 
is  obtained  by  considering  the  action  of  convective  fluid  mo- 
tions as  they  deform  the  toroidal  field.   The  convective  process 
is  pictured  to  be  an  arbitrary  distribution  of  ra-dially  directed 
rising  and  sinlcing  columns  of  fluid.   The  earth's  rotation  en- 
ters into  the  process  a  second  time  v'hen  the  Coriolis  force 
imparts  a  cyclonic  motion  to  each  convective  column.  For  the 
case  of  a.n  cutwa.rdly  directed  column  passing  through  a  small 
region  of  the  toroidal  field,  each  line  of  force  is  lifted  and 
then  twisted  90°  counterclockv/ise  to  form  an  elemental  loop  of  ' 
magnetic  flux  lying  in  a  meridional  plane. (Fig.  6).  Parker 
further  showed,  that  for  any  arbitrary  distribution  of  convec- 
tive columns,  the  elemental  loops  of  flux  will  coalesce  and 
regenerate  the  poloidal  field. 
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FIG.  5.     GENERATION     BY     NON-UNIFORM      ROTATION     OF 
A     TOROIDAL     FIELD     FROM    A    POLOIDAL    FIELD. 
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FIG.    €.     GENERATION    OF  Nr^ERlDlGMAL    LOOPS    Or     ^iAGNETIC     FLUX 
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THSCRI2:S  C?  THE  EARTH'S  EXTERHAL  FIELD 
Tiieorj  of  the  Solar  Dail^''  Variation 

As  the  prinary  source  of  tlie  S  or  solar  dally  variation, 
Stewart  proposed,  "that  convective  currents  established  by  the 
sun' s  heating  influence  in  the  upper  regions  of  the  atmosphere 
are  to  be  renarded  s.s   conductors  moving  across  lines  of  mag- 
netic force,  and  thus  the  veliicle  of  electric  currents  which 
act  upon  the  magnet".   This  atmospheric  djTiamo  theory  for  the  S 
variation  was  developed  into  a  c_uantitative  form  by  Schuster 
(1908),  and  reriains  quite  well  accepted  tod.a3^.   Other  theories 
involving  atriosphcric  thermoelectric  currents  or  the  change  of 
paramagnetic  behavior  of  oxygen  when  heated  do  not  have  any 
basis  in  ph3^sical  fact. 

The  djT.amo  action  given  by  Schuster  tahes  place  in  a  thin 
spherical  shell  of  the  ionosphere.   The  d^Tnamo  field  is  given 
by 

E  =  V  X  K  (er.u  system)  (44) 

■where  v  is  the  velocity  vector  in  the  T)resence  of  a  magnetic 
field  K.   If  the  r:Otion  of  the  air  of  the  shell  is  supposed 
horizontal  and  irrotational,  it  is  associs-ted  with  a  velocity 
potential 'y,  such  that, 

V  zr  -  gradT  ,  (45)' 

Let  the  x  direction  be  southward,  the  y  direction  be  eastward, 
ejnc.   the  z  direction  be  upward,  and  let  e.   be  the  radius  of  the 
shell.   Since  v^  =  0,  "Vis  independent  of  z,  and  equations  (44) 
and  (45)  can  be  combined  to  give 
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E 


X 


Ey 


-H. 


sin  e  ^0 


Hz  ^^^ 


a  ^?e 


(46) 


v/here  e  and  0   arc  the  angles  of  co-latitude  and  of  east  lon£;i- 
tude.   Thus  the  d^mamo  field  depends  onljr  on  the  vertical 
component  cf  the  G^o^is-S^^tic  field. 

The  existence  of  a  steady"-  current  system,  v/ithcut  a  c;roi;- 
ing  accumulation  of  charge,  implies  that  the  scuthv;ard  and 
eastv.'ard  components  (U,  V)  of  the  current  densit3''  are  associ- 
ated v/ith  a  current  function  R.   Therefore, 


U  -  hR 

a  sin  e  ^0 

V  =  -  ^R 


(^7) 


bQ 


The  electric  force  required  to  drive  thds  current  must  have  the 
components  U/iC,  and  V/K,  v.'here  K  is  the  conductivitj'"  of  the 
spherical  shell.   Kov:ever,  the  dynamo  components  Ex  and  Ey  are 
not,  in  general,  related  to  any  function  R  in  the  manner  of 
U/K  and  V/K.   The  dyna,mo,  field  will  set  up  a  charge  distribu- 
tion whose  field,  together  with  the  d3rnamo  field  and  an}'-  addi- 
tional field  due  to  self  or  mutual  induction,  is  such  that  a 
current  system  having  the  current  function  R  is  produced.   If 
self  or  mutual  induction  is  ignored,  the  equation  satisfying  R 
is 


^^^    +  A  (sin  e  Q^] 
sin  e  h0'^       <je  '      ^Q  . 


^'^    ^H    +  sin  e  ^  aR 
J^   sin  e  >)jZ)        ^e  W 


-    B-Y.^ 


<^(vxKz)  +  MvyHz  sin  9) 


(48) 


Thus  R  is  determined  if  Hz,  K,  and  v  (or  the  velocity  potential 


A^3 


^)  are  l-inovni. 

TIiecrie:3  cj?  I-la^nstic  Storns 

Aii5''  theory  of  nsLGnetic  Gtonv.s  must  c::plain  not  cni3''  the 
masnetlc  disturbance  variation  D,  "but  also  such  allied  phenom- 
ena as  the  auroras  and  the  Van  Allen  "belts.   These  phenomena 
i-rill  be  briefly  considered  before  the  discussion  of  magnetic 
storm  theories.  Auroral  displays  occur  in  t\io   primary  zones 
20°  to  25°  from  the  geomagnetic  poles;  Vestine  (I960)  has  sho:-ni 
that  these  zones  are  linlied  hj   geomagnetic  linss  of  force.   The 
displays  alwa^^-s  ts.ke  place  at  timios  \ihen   the  m.agnctic  field  is 
disturbed,  but  not  all  storms  are  accompanied  by  auroras.  Dur- 
ing m-ore  intense  disturbances,  the  auroral  zones  tend  to  move 
tovjard  Icvrer  latitudes.   I-Ioct  of  the  displaj^s  are  at  a  height 
of  approximately  100  hm..   Since  the  spectra-  of  the  displays 
have  many  strong  hydrogen  lines,  proton  bomba2?dm.ent  of  the  upper 
atm:osphere  must  certainly  play  a  part  in  aurora  generation.   The 
incidence  of  auroras  is  at  a  maximum  one  or  ti70  years  after  the 
•peah  of  the  eleven  year  sunspct  c^^'cle.   For  any  year  auroras  are 
miore  prevalent  during  the  equinoctial  season. 

Regions  of  intehse  radiation  in  the  space  surrounding  the 
earth  were  first  discovered  oj   Van  Allen  et  al.  (195-^)  through 
analysis  of  data  obtained  by  earth  satellites  Explorers  I  and 
III.   This  data  and  the  data  from  later  satellites  showed  the 
existence  of  t-'.jo  distinct,  widely  separated,  high  intensity 
zones  of  trapped  particles 'encircling  the  earth.  T.ie   trapped 
radiation  is  largely  protons  and  electrons  that  are  constrained 
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uo  uroclioidal  orbits  by  the  geomagnetic  field;  tlie  p^.rticles 
ovePitually  leal:  out  of  tlie  radiation  zones  by  virtue  of  cclli- 
cional  enerG2^  loss  and  scattering  in  tlie  thin  outer  atmosphere. 
The  inner  Van  Allen  belt  lies  1.5  earth  radii  from  the  ceo- 
center;  its  particles  seen  to  be  the  products  of  radioactive 
decay  of  neutrons,  these  neutrons  beinj;;  released  fron  cosnic 
ra.y  induced  nuclear  disintegrations  in  the  top  of  the  atmos- 
phere.  The  cciter  Van  Allen  belt  is  centered  at  3  to  4  earth 
radii  from  the  ceocenter;  the  source  of  its  particles  appears 
to  be  solar  plasma  that  has  reaciied  the  earth.   This  outer  belt 
may  play  a  central  role  in  auroras,  geom^agnetic  storms,  and  at- 
mospheric heating.  An  artificial  Van  Allen  belt  iras  created 
by  the  A-rgiis  experiment  in  195>-';  three  small  nuclear  devices 
were  detonated  at  high  altitudes  to  inject  electrons  into 
trapped  orbits  in  the  earth's  mxagnetic  field  (Christof ilos, 

Birkels-nd  (1896)  gave  the  first  promising  suggestion  on 
the  correlation  of  auroras  and  magnetic  disturbances.   Ke  sh.owed 
that  I'hen  an  electron  stream,  was  projected  towards  a  magnetized 
sphere  in  an  evacuated  chamber,  the  electrons  were  deflected 
towards  the  poles  oj   the  m:agnetic  field.  A  behavior  quite  sim.- 
ilar  to  the  earth's  aurora  took  place;  two  circule.r  discharge 
zones,  one  about  each  pole,  appeared  upon  the  sphere. 

Storm.er  (19C7,  1911,  1912)  developed  m.athematically  the 
results  of  Birkeland's  experiments  by  calculating  the  possible 
paths  of  particles  of  one  sign  m.oving  under  the  influence  of 
the  geomxagnetic  field,  ^fith  the  assum-ption  that  these  particles 
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vrere  emitted  fro~:  the  sun,  Stcrmcr  ires   able  to  sliov:  tha,t  zuch 
particles  coLild  reach  the  earth'  s  upper  atmosphere  only  in  the 
auroral  zones.  A  further  consequence  of  these  calcul2.tions 
i/as  the  prediction  of  a  ^Testward  flo^.Ting  ring  current  about 
tlie  earth.   The  calculations  showed  that,  in  a 'bundle  of  solar 
particles  reaching  the  earth,  many  particles  would  be  deflected 
around  the  earth  in  the  equatorial  plane.   The  resultant  ring 
current  could  account  for  the  main  phase  decrease  in  K  during 
a  mar;;netic  storm  and  could  also  explain  the  tendency  of  the 
auroral  zones  to  shJ.ft  equatorwards  during  a   severe  magnetic 
disturbance. 

One  fatal  objection  to  the  Birlieland-Storner  theory  is 
that  if  the  emitted  corpuscles  from  the  sun  were  of  one  sign 
only,  they  would  disperse  through  mutual  electrostatic  repul- 
sion long  before  reaching  the  earth.   Nevertheless,  the  Stcrmer 
motion  of  charged  particles  pla3i's  a  predominant  role  in  the 
geners-tion  of  the  auroras  and  the  Van  Allen  belts,  , 

Cha.pman  and  Ferraro  (1931)  hypothesized  that  the  effects 
of  a  magnetic  storm  were  due  to  an  electrically  neutral  stream 
of  sola.r  particles.   The  details  of  a  sola.r  process  that  could 
produce  such  a  neutral  ionized  stream  are  not  yet  completely 
knovfxi.   The  time  interval  of  about  a  da^r  between  the  occurrence 
of  a  solar  flare  and  the  onset  of  a  mzagnetic  storm  requires 
that  the  stream  have  a  velocity  of  6C0  to  14CC  miles  per  sec- 
ond.  It  vras  estimated  that  the  width  of  the  stream  front  was 
at  least  lo'^  hm.,  so  that  the  surface  of  tlie  front  appears  as 
a  plane  as  it  nears  the  earth.   The  comDosition  of  the  stream. 
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vjas  considered  to  be  larcel^^  electrons  and  protons,  irlfh   ne^li- 
[ylole   recombination  taking  place  during  tlie  passable  of  the 
stream  from  tlie  sun  to  tlie  earth. 

As  this  plane  front  of  ionized,  hi^hl^;^  conducting  gs'-s  ap- 
proaches the  earth,  currents  are  induced  in  the  stream's  surface 
by  the  QecmaG'^otic  field.   The  resultant  magnetic  effect  is 
equivalent  to  e/r.   image  dipole  situated  inside  the  stream  at  the 
S32iie  distance  from  the  stream  surface  as  the  earth  dipole  is 
outside  the  stream.   Such  an  image  dipole  compresses  the  lines 
of  force  of  the  earth  dipole  on  the  side  nearest  the  stream., 
a,nd  causes  an  increase  in  r;agnetic  intensity  typical  of  the 
first  "ohase  of  a  magnetic  storm. 

In  svreeping  by  the  earth,  the  solar  stream  is  repelled 
from  the  earth's  surface  b^r  the  compressed  geomagnetic  field 
lines,  and  a  hollov;  space  forms  in  the  stres.m  (Fig.  7).   Sur- 
face charges  on  the  vralls  of  the  holloi;  produce  an  electric 
field  that  gives  rise  to  a  discharge  across  the  intervening 
gap.   This  discharge  leads  to  a  net  i'estir3.rd  flo'-r  of  current 
around  the  earth.  \Ilien   the  stream  has  passed  bej-'ond  the  earth, 
a  ring  current  ••.'•ill  be  left  behind.   This  current  produces  the 
decrease  in  K  charade teris tic  of  the  m.ain  phase  of  a  magnetic 
storm. 

Alfv6n  (1955)  accepted  the  idea  that  the  solar  stream  con- 
sists of  a  neutre^l  ionized  gas,  but  he  proposed  also  that  the 
stream,  because  of  its  high  electrical  conductivity,  may  pos- 
sess e.ii   embedded  magnetic  field.   Thus,  as  the  stream  moves  ''oy 
a  point,  an  electrical  field,  E  =  v  X-H,  ■'.•.all  appear;  in  effect 
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tlie  strean  is  electrical  13?-  pola.rized.   By  aeohods  developed  "bj 
Storr.er,  Alfven  calculs^ted  tl-.e  nction  of  electrons  under-  the 
influence  of  this  electric  field  and  the  earth  s  dipole  field. 
These  electrons  drift  in  trcchcidal  orbits  about  the  earth;  the 
motion  is  equivalent  to  a  rino;  current  in  the  eG_uatcrial  ^jlane. 
D"arins  the  onset  of  a  magnetic  storm  Alf v^n'  s  rins  current 
flows  west  to  east;  the  current  reverses  in  direction  during 
the  main  phase  of  the  storm.   In  ccmparins  his  ring  current 
with  the  rin5  current  of  the  Chapman-Ferraro  theory,  Alfv(=;n 
pointed  cut  that  an  electron  moving  in  a  trochoidal  path  which 
has  an  overall  circular  drift  pattern  represents  a  stable  orbit, 
whereas  it  had  never  been  proved  that  a  Ginple  ring  current  was 
stable. 

Still  another  theor3^  c>f  the  magnetic  storm  V7as  developed 
by  Singer  (1957).   This  theory  considers  that  the  stream  from 
a  solar  flare  gives  soriO  of  its  energy  to  the  production  of  a 
shoch  wave  in  the  interplanetary  gas.   This  shocli  v.'ave  propo- 
gates  at  a  higher  velocit^^  than  the  solar  flare  itself,  and 
a.rrives  at  the  earth  up  to  nine  hours  earlier  than  the  solar 
stream,   ''vhen  the  shcch  wave  reaches  the  earth,  it  is  retarded 
in  the  equatorial  plane  by  the  geomagnetic  field.  Above  and 
below  tills  plane  the  shcch  wa.ve  is  bent  and  gLiided  "oj   the  lines 
of  force,  and  becomes  compressed  'b'j   the  field  as  it  travels  into 
the  atmosphere  of  the  auroral  zone.   There  auroral  displays  are 
produced  and  a  current  is  initiated  which  is  responsible  for 
the  sudden  commencem.ent  of  a  magnetic  storm;.. 

The  main  -ch.ase  of  the  magnetic  storm  results  when  the 
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solai''  strcar:  particles  reach  the  ec.rtn.   The  solar  £:as  pene- 
trates into  the  ■::eoiiiasnetic  field,  which  has  been  perturbed  by 
the  shock  ija-VG.  After  beinG  scattered,  the  charged  particles 
are  attaclied  tc  lines  cf  na^rnetic  force  and  exGciite  spirc-lin^ 
or  trcchoidal  notions  a/ocait  tliese  lines.  At  the  sane  tiiie 
these  trapped  particles  drift  in  azimuth  around  the  earth  and 
generate  a  large  uest'-ard  flooring  current  of  the  right  magnitude 
tc  cause  the  iiain  phase  effect  of  the  stern. 

Singer's  concept  cf  a   trapped  region  of  solar  particles 
seems  amply  supported  by  the  discovery  of  the  Van  Allen  belts, 
and  by  the  results  of  Project  Argus.   The  ps-rticles  oscillate 
rapidly  to  and  fro  betvjeen  rL£.gnetic  mirror  points  at  fairly 
high  northern  and  southern  1  altitudes.   Simultaneously,  the^r 
circle  afoout  the  ms.gnetic  field  lines,  and  also  drift  arcund 
the  earth,  protons  west^-jard,  and  electrons  eastward. 

Another  recent  ma.gnetic  storm  theory  was  presented  hj 
Dessler  and   Parker  (1959).   Iii  this  theory  the  storm  is  caused 
by  the  collision  of  the  earth'  s  m.agnetic  field  i;ith  relatively 
dense  plasma  clouds  ejected  from,  the  sun.   The  initial  phase 
of  the  storm  occurs  -vihen   a  solar  plasma  cloud  pushes  s.gainst 
the  earth's  field  and  compresses  it.   The  rapid  increase  in 
magnetic  intensity  at  the  beginning  of  the  storm  is  propogated 
rapidljr  around  the  earth  by  a  hydromagnetic  wave. 

Upon  diffusion  of  the  pla,sm.a  into  the  geomagnetic  field, 
the  main  phase  begins.   The  trapped  particles,  mostly  protons, 
exert  tliree  stresses  upon  the  geomagnetic  field.   The  chief 
stress,  due  tc  centrifugal  force  of  the  particles,  accounts 


for  most  of  the  nain  phase  changes.   The  recovery  frcn  the  main 
phase  is  attributed  to  the  relief  of  the  geomagnetic  field 
stress  by  a  transfer  of  the  energy  of  the  trapped  protons  to 
atoms  of  neutral  hydrogen  thrcugli  an  ion-atom  cha^rge  exchange. 

CCICCLUSION 

At  the  earth's  surface  the  geomagnetic  field  is  derivable 
from  a  scalar  potential.   3y  the  use  of  spherical  harmonic 
analysis  it  is  possible  to  determine  the  components  of  the 
field  due  to  external  and  internal  sources.   The  internal 
field,  much  larger  than  the  external  field,  has  a  large  compo- 
nent that  can  be  produced  'oy   a  simple  dipole,  and  a  smaller 
non-dipole  or  residual  component.   With  respect  to  the  earth's 
dipole,  a  system  of  geomagnetic  co-ordinates  useful  in  the 
study  of  magnetic  storms  can  be  derived. 

The  eo.rth's  magnetism  exhibits  tivo  distinct  tj^pes  of  var- 
iation, secular  and  transient.   The  secular  variation,  of  in- 
terns-l  origin,  is  clearly  associated  v:ith  the  non-dipole  field, 
and  has  displayed  s.   vjestv.-ard  drift  during  the  past  half  century, 
Three  major  forms  of  transient  varis.tion  exist:  (1)  the  solar 
daily  variation,  (2)  the  lunar  daily  variation,  and  (5)  the 
magnetic  disturbance  variation.  All  thjr^ee  form.s  are  largely 
of  external  origin. 

For  a  sim.ple  m-odel  of  the  earth's  miagnetism,  the  internal 
field  dipole  m&j   be  placed  at  the  center  of  the  earth.  A  sorie- 
what  better  representation  is  possible  by  displacing  the  dipole 
25c  km.  from  the  geocenter.   The  residual  field  r.ay  be  repre- 
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sented  "by  a  distribution  of  dipoles  or  current  functions  at 
the  surface  of  the  earth's  liquid  core.   A  combination  of  an 
overhead  ionospheric  cur-rent  system  and  an  equatorial  ring 
current  produces  a.n  adequate  model  of  the  transient  magnetic 
variations.  An  accurate  simulation  of  the  geomagnetic  field 
is  obtained  by  combining  the  models  for  the  dipole  field,  the 
residual  field,  and  the  external  field. 

A  dynamo  process  seems  the  most  likely  source  of  the  in- 
ternal field.   The  non-linearity  of  the  basic  magne to-hydro - 
dynam.ic  equations  allov:s  onl3'-  a  partial  solution  of  the  dynamo 
action.   For  certain  assumed  fluid  motions,  a,  solution  satis- 
fying Maxv/ell's  oquations  and  the  postulated  boundary  conditions 
does  exist,  but  the  solution  is  not  necessarily?"  unique.   This 
solution  requires  an  internal  toroidal  field  that  is  not  de- 
tectable at  the  earth' s  surface  and  is  much  stronger  than  the • 
earth's  dipole  field. 

Tides  in  the  upper  atmosphere  appear  responsible  for  the 
solar  d3.ily  variation  in  the  m.agnetic  field.  Emitted  solar 
particles  reaching  the  earth  lead  to  m.agnetic  storms  and 
8.uroral  displays.   These  p8,rticles  also  enter  into  the  form.a- 
tion  of  the  Ve-n  Al-  i  belts.   The  development  of  a  comprehensive 
magnetic  storm  theory  must  certainly  av/ait  a  more  complete 
measurement  of  the  external  geom.agnetic  field  and  the  solar 
i;ind  of  em.itted  plasmia  particles. 
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This  report  ^ivec  a  detailod  descripticn  of  tl-.e  ^eoir.as- 
notic  field,  presents  some  r-^hyslcal  iQcdels  of  the  field,  and 
discusses  a  numter  of  theories  relating  to  the  field' s  origin. 

The  magnetic  elements  used  in  geomagnetic  surface  measure- 
ments are  defined,   With  the  use  of  spherical  hamonic  analysis 
a  scalar  magnetic  ^potential  function  naj/  "be  obtained  fron  \;orld 
i;ide  magnetic  data.   This  potential  function  is  separable  into 
tv;o  component  potentials,  one  due  to  internal  sources,  and  the 
other  due  to  externa-1  sources. 

The, internal  sources  produce  a  much  larger  field  than  the 
external  sources;  this  internal  field  has  e,   large  component 
that  could  be  produced  "oj   a  sim.ple  dipole,  and  a,  smaller  non- 
dipole  or  residual  component.   The  residual  internal  field  is 
associated  Ytxth   the  secular  variation  of  the  earth's  field. 

External  sources  superimpose  upon  the  earth's  mean  field 
three  major  forms  of  transient  m.agnetic  variation:  the  solar 
daily  variation,  the  lunar  dail^y  vs-riaticn,  and  the  magnetic 
disturbance  variation. 

The  simplest  model  of  the  earth's  m.agnetism  V70uld  repre- 
sent the  main  field  bj^  either  a  centered  dipole  or  "oj   a  slightly 
eccentric  dipole.   hore  detailed  models  must  also  include  the 
internal  non-dipole  field  and  the  externa.l  field.   A  distribu- 
tion of  dipoles  or  current  functions  at  the  surface  of  the 
earth's  core  could  generate  the  non-dipole  field.  A  combination 
of  an  ionospheric  current  s^'-stem  and  an.  external  ring  current  in 
the  equatorial  plane  produces'  an  adequate  m.odel  of  the  external 
field. 


Various  theories  on  the  cri^ij^  c^  t^^-S'  internal  field, 
oaced  en  ferrcLiagnetisn,  spontaneous  magnetization  of  rotating 
bodies,  syste:r.s  of  rcta.tinr;  electric  charge,  or  slov.'ly  decaying 
current  systens,  are  shov.'n  to  fail  in  nany  aspects.   The  in- 
ternal field  is  probabl3''  due  to  internal  current  s^^stems  main- 
tained either  by  the  thermoelectric  effect  or  by  a  dynamo 
action;  the  latter  process  seens  more  likely,  as  it  requires 
fe-.fer  assurrptions. 

The  basic  equations  governing  d3rnai:o  action  in  the  earth'  s 
fluid  core  are  non-linear  and,  therefore,  not  easily  amenable 
to  analT'sis..  Some  qualitative  results  m.a^''  be  obtained  by  a 
dimensional  analysis  of  those  equations.  Another  approach  to 
the  problem:  c-ssumed  a  particular  fluid  motion,  and  then  showed 
th£it  a  self  regenerating  d3rLiaiiio  action,  consistent  ■-.'•ith  x-Iazvell' s 
equations,  is  possible. 

The  solar  daily  variation  in  the  magnetic  field  is  m.ost 
lil-iely  a  result  of  upper  atriospheric  tides.   The  rxagnetic  dis- 
turbance varla-tion  and  its  associated  phenomena  apparently  re- 
sult from  the  interaction,  of  emitted  solar  particles  uith  the 
earth's  field.  Any  comprehensive  magnetic  storm  theory  m.ust 
airait  a  m.ore  com.plete  determination  of  the  external  geomagnetic 
field  and  the  nature  of  the  solar  v;ind. 


